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I. INTRODUCTION 



Generally, in coastal and ocean waters, the velocity profiles are typically established by 
bottom friction and by surface wind stress and so are varying with depth. Currents generate 
shear at the bed of the sea or of a river. For example ebb and flood currents due to the tide 
may have an important effect on waves and wave packets. In any region where the wind is 
blowing there is a surface drift of the water and water waves are particularly sensitive to 
the velocity in the surface layer. 

Surface water waves propagating steadily on a rotational current have been studied by 
many authors. Among them, one can cite Tsaoi, Dalrymple 2 , Brevik-, Simmen & Safmann^, 
Teles da Silva & Peregrine^, Kishida & Sobey-, Pak & Chow 7 , Constantin^, etc. For a 
general description of the problem of waves on current, the reader is referred to reviews 
by Peregrine 9 , Jonsson 10 and Thomas & KlopmaiA On the contrary, the modulational 
instability or the Benjamin- Feir instability of progressive waves in the presence of vorticity 
has been poorly investigated. Using the method of multiple scales Johnson^ examined the 
slow modulation of a harmonic wave moving over the surface of a two dimensional flow 
of arbitary vorticity. He derived a nonlinear Schrodinger equation (NLS equation) with 
coefficients that depend, in a complicated way, on the shear and gave the condition of linear 
stability of the nonlinear plane wave solution by writing that the product of the dispersive 
and nonlinear coefficients of the NLS equation is negative. He did not develop a detailed 
stability analysis as a function of the vorticity and depth. Oikawa, Chow & Benney^ 
considered the instability properties of weakly nonlinear wave packets to three dimensional 
disturbances in the presence of shear. Their system of equations reduces to the familiar 
NLS equation when confining the evolution to be purely two dimensional. They illustrated 
their stability analysis for the case of a linear shear. Within the framework of deep water 
Li, Hui & DonelaniA studied the side-band instability of a Stokes wave train in uniform 
velocity shear. The coefficient of the nonlinear term of the NLS equation they derived was 
erroneous as noted by Baumstein^. The latter author investigated the effect of piecewise- 
linear velocity profiles in water of infinite depth on side-band instability of a finite-amplitude 
gravity wave. The coefficients of the NLS equation he derived were computed numerically 
because he did not give their expression as a function of the vorticity and depth of the shear 
layer, explicitly. Instead, he calculated these coefficients for specific values of the vorticity 
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and depth of shear layer. Choi^ considered the Benjamin- Feir instability of a modulated 
wave train in both positive and negative shear currents within the framework of the fully 
nonlinear water wave equations. For a fixed wave steepness, he compared his results with 
the irrotational case and found that the envelope of the modulated wave train grows faster 
in a positive shear current and slower in a negative shear current. Using the fully nonlinear 
equations, Okamura & Oikawa^ investigated numerically some instability characteristics of 
two-dimensional finite amplitude surface waves on a linear shearing flow to three-dimensional 
infinitesimal rotational disturbances. 

The present study deals with the modulational instability of one dimensional, periodic 
water waves propagating on a vertically uniform shear current. We assume that the shear 
current has been produced by external effects and that the fluid is inviscid. In section [TT] 
a NLS equation (vor-NLS equation) for surface waves propagating on finite depth in the 
presence of non zero constant vorticity is derived by using the method of multiple scales. 
In subsection III CI it is shown that the heuristic method to derive a NLS equation from a 
nonlinear dispersion relation is not valid when vorticity is present. This is a consequence of 
the coupling between the mean flow due to the modulation and the vorticity. Section IIHI 
is devoted to a detailed stability analysis of a weakly nonlinear wave train as a function 
of the parameter kh where k is the carrier wavenumber and h the depth and of vorticity 
magnitude. Consequences on the Benjamin- Feir index are considered, too and a conclusion 
is given in section IIVI 

II. DERIVATION OF THE VOR-NLS EQUATION 

The undisturbed flow is a weakly nonlinear Stokes wave train propagating steadily on a 
shear current that varies linearly in the vertical direction y. The wave train moves along the 
x-axis. The y-axis is oriented upward, and gravity downward. Naturally, i and j are unit 
vectors along Ox and Oy. When computing vector products, we shall also use k = i A j. 
The depth h is constant and the bed is located at y = — h. Let Q be the magnitude of the 
shear. There is a potential ip(x,y,t) such that the velocity writes 

V = Qyi + V(p(x,y,t) (1) 

since for a two dimensional flow of an inviscid and incompressible fluid with external forces 
deriving from a potential the Kelvin theorem states that the vorticity is conserved. The 
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Figure 1. Shear flow with Q > Figure 2. Shear flow with f2 < 

(waves propagating downstream) (waves propagating upstream) 

variable t is the time and — Q is the vorticity in all the fluid that can be negative or positive 
as illustrated in figures [I] and [2j respectively. Note that the reference frame is in uniform 
translation with regard to that of the laboratory. Hence, the velocity of the undisturbed 
flow vanishes at the surface. 



A. Governing equations 

As the perturbation is assumed potential, the incompressibility condition V.V^ = implies 
that the velocity potential satisfies the Laplace's equation 

Aip = - h< y < r)(x,t) (2) 

The fluid is inviscid and so the Euler's equation writes : 

V{^> t + \Y 2 + - + gy) = VAlu (3) 
z p 

where u is the vorticity vector : u = —Qh. 

We introduce the stream function ip associated to the velocity potential through the Cauchy- 
Riemann relations : 

4>y = fx, ^x = (4) 

We notice that 

VAw = V(^V + ^) (5) 
so, we can rewrite equation (j3J) as follows 

Vfat + \& + + Vy<t> x + - + gy - = 0. (6) 
2 2 y p 

4 



This equation may be integrated once : 

IIP 

ft + -jp\ + o¥>l + My<t>x + - + gy-ni(> = f(t) (7) 
2 2 y p 

We write this equation at the free surface of the fluid. The notation $ means that <p is 

calculated on the free surface. The same convention is used for the derivatives of </? or if). 

This convention will be used in the whole paper. 

The pressure on the free surface is the atmospheric pressure that can be considered as a 
constant, and incorporated in the RHS of equation (J7J). 

It is possible to add to the velocity potential function a primitive of the right hand 
side f(t) of this equation, so that this term vanishes. The equation becomes 

$ t + U>1 + l -& y + n^ x + gv -nq = o. (8) 

The kinematic condition is written as follows 

% = Vx(®x + ton) + Vt (9) 



The governing equations are then 

VV = 0, -h<y<r){x,t) (10) 

<p y = o, y = -h (11) 

th + ($x + torfiv* ~ % = (12) 

$t + W x + l$l + nr]$ x + gTi-m = (13) 



To reduce the number of dependent variables, we derive the dynamic condition with respect 
to x and we use the Cauchy-Riemann conditions to eliminate the stream function. The 
result is 

'I'o- + ®tyVx + $x{$xx + ®xyVx) + ®y($xy + ®yyVx) 

(14) 

+nrj x <l> x + nrj($ xx + $ xy r} x ) + gr] x + fi($ y - <$> x -q x ) = 
Equations (fT0"j) - ([T3"j) are invariant under the following transformations : if — > —if, t — > —t, 
Q — > —Q and ^ — > — Hence, there is no loss of generality if the study is restricted 
to waves with positive phase speeds so long as both positive and negative values of Q are 
considered. 
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B. The multiple scale analysis 

We seek an asymptotic solution in the following form 

+00 +00 
if = (p n exp[in(kx — cot)], r/ = r) n exp[in(kx — ut)] (15) 

7i=— 00 n=— 00 

where k is the wavenumber of the carrier and u its frequency. 

We assume = ip* n and r\_ n = 77* so that tp and r\ are real functions. 

Then ip n and r\ n are written in perturbation series 

+00 +00 

¥n = eVnj, = e^ry (16) 

where the small parameter e is the wave steepness. 
We assume ipoo = and r] 00 = 0. 

Following Davey & Stewartson^, we consider a solution that is modulated on the slow time 
scale r = e 2 t and slow space scale £ = e(x — c g t), where c g is the group velocity of the carrier 
wave. 

The new system of governing equations is 



eVtt + <Pw = °> ~h < y < rj(C, T ) (17) 

^ = 0, y = (18) 

e 2 ?7 r - ec 9 ?7 f + e 2 $ ? ?7 ? + eft^ - $ y = (19) 

+ + e%^y + eQyQyyriz + (Hh/.A'c + ( 2 < > < I > , i (20) 
+ e 2 ^?7$ ?y ?7 € + e#?7 C + fi$„ - e 2 fi$ € r^ = 
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Substituting the expansions for the potential ip into the Laplace equation and using the 
method of multiple scales we obtain 

Voiyy = (21) 

-£;Vll + IfiUyy = (22) 

(f02yy = (23) 

-k 2 tfi 12 + 2ik(pu£ + IfiUyy = (24) 

-4fc V22 + V^2 2ra = (25) 

Voice + V?03 TO = ( 26 ) 

-ArVi3 + 2?A;v9i2 5 + + V?i3 TO = (27) 

-4&V23 + 4ifcv?22C + V- , 23j/j/ = (28) 

-9A;V33 + V?33 TO = (29) 

Solving these equations and considering the bottom conditions we obtain : 

VOly = (30) 

V02y = (31) 

cosh[fc(y + h)} 

^ = A cosh(kh) (32) 
_ D cosh[k{y + h)} {y + h) sinh[fc(y + h)} - ha cosh[fc(y + h)} 
12 cosh(A;/i) cosh(fc/i) 
cosh[2fe(y + ZQ] 

^ 22 " F cosh(2^) (34) 

Vo3j, = -(?/ + M0oi& ( 35 ) 

(y + /i) sinh[/c(y + /i)] 



V?i3 = (/icrA^ - i^)- 



cosh(/c/i) 



^- 2 /-, ~ ,,,,,,, ( ., „ . (v + /i) 2 N cosh[Hy + h)] 
+ (B + —(l-2tanh 2 (kh))Att + thaD{:-Att- 



(36) 



2 v v " K ? ?? 2 y cosh(A;/i) 

The next tedious step is to use the relations obtained from the kinematic and dynamic 
conditions. Let us set Vl = ^. Herein, it is important to emphasize that this parameter 
does not correspond to a dimensionless vorticity because the frequency u depends on Q. 
Furthermore, we set X = crQ, where a = tanh(/c/i), because this term will occur many times 
in the following polynomial expressions. 

• Terms in eE° : They give no supplementary information 



Terms in eE 1 : They give a linear dispersion relation 



kc 2 p + a(c p n -g) = (37) 

From this linear dispersion relation it is easy to demonstrate that we have always 
X > -1 or H > -1/(7. 

We also get a relation between the fundamental modes of the velocity potential at the 
surface and free surface elevation : 

mi = i-A=—(l + X)A (38) 
Cp 9 

where c p = u/k. 

From the linear dispersion relation the group velocity is : 

c (l-a 2 )kh + a(l + X) 
Ca ~a X 2 + X ^ 

We recall that X > — 1, so there is no singularity in this expression. 
Terms in e 2 E° : They only give, after simplifications : 

Von = (40) 

Terms in e 2 E 1 : They give a system of two equations with two indeterminate coeffi- 
cients i] i and i]i2 that can be found after some calculations : 

Voi = (41) 

and 

Vi2 = - g \°9 + Ml " <J 2 MAz + j (1 + X)D (42) 

Terms in e 2 E 2 : They give a system of two linear equations with F and r] 2 2 as un- 
knowns : 

F^ ^-^^ A' (43) 



and 



V22 = -^-[3 - a 2 + (3 + a 2 )X + X 2 ]A 2 (44) 
2c p a 
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Terms in e E : The first-order mean flow can be obtained from the following expres- 
sion 

1 2 



[c g (c g + tth) - gh]ip m 
and 



9 -^{2 + X) + k 2 c g {l-a 2 ) 



\A\< (45) 



9V02 = (c g + nh)<P m - k\l - a 2 )\A\ z (46) 

Terms in e 3 E 1 : We derive two equations from which it is possible, after tedious 
computations, to eliminate 7713. The coefficients B and vanish owing to the linear 
dispersion relation. The remaining terms are : A time derivative, a dispersive term, 
a nonlinear term and a term involving the mean flow that we can substitute by its 
expression taken from the other equation. Finally a nonlinear Schrodinger equation 
with vorticity is derived (the vor-NLS equation) 



where 



iA T + LA# = P I A | 2 A (47) 



L = k 2 a{ 2 + X f {1 ~ a2)[1 ~ ^ + (1 " P)X] ~ ° p2 (48) 

P - fc % ( U + VW)= k ^ U + VW ) (49) 
2(2 + X)ga* \ * ) 2 (1 + X) (2 + X)ua* 1 J 



3 (l + X)(2 + X)+p(l-^ 



with 



(50) 



U = 9- 12a 2 + 13a 4 - 2a 6 + (27 - 18a 2 + 15a 4 )X 

+ (33 - 3a 2 + 4a 4 )X 2 + (21 + 5a 2 )X 3 + (7 + 2a 2 )X 4 + X 5 
V = {l + Xf[l + p + fiU) + 1 + X - pa 2 - paX (51) 



W = 2a 6 — /v _ ' rK - (52) 

ap(p + pQ) - p{l + X) 



p = kh (53) 
a = tanh(//) (54) 



p = — (not to be confused with the density) (55) 

Cp 

The relation ( I38p permits to replace the velocity potential A by the elevation a. 

ia T + La,££ = M I a \ 2 a (56) 
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where a is the envelope of the surface elevation and 

M = 8(1 + X)(2 + XK (5?) 
C. The case of infinite depth 

Let us discuss what happens when depth goes to infinity in order to compare our results 
to those of Li et al. 14 or Baumstein 1 ^. Moreover, we shall show the importance of the 
coupling between the mean flow and vorticity at third order. At this order before deriving 
equation (j47p the following coupled equations empasize the coupling between the mean 
flow O i£ an d vorticity Q. 

k 3 cl 



-p 



9° 



(1 + of]) 2 (c p + c g + khc p Q) + c p (l + o-O) 
— (c g + c p khQ)a 2 ~\ 4>q\^A — iu(2 + aQ)A T 



+ {c 2 -gh + gha [a + kh(l - o 2 )] + c p kh(l - o 2 ) (c g - c p kha) 0} A. 

k 5 c 2 r — 
" 9 - 12o 2 + 13o 4 - 2o 6 + (27 - 18o 2 + 15o 4 )ofi 



(58) 



+ (33 - 3o 2 + 4o 4 )o 2 2 + (21 + 5o 2 )o 3 3 



+ (7 + 2o 2 )o 4 4 + a 5 Q 5 



\A\ 2 A 



with 



co (2 + oQ) + k 2 c p Cg(l — o 



2 

2 



M (59) 



° 1? c p [c g (c g + Qh) - gh] 

The mean flow 0oi£ verifies 

gk(2 + Tl) , |2 

Jim h^ = ^- L A 2 60 

Cp(\lCg — g) 

The coefficient of the mean flow, induced by the modulation of the envelope, in equa- 
tion (1581) is of order 0(h) so that the product has a finite limit when h — > 00. More precisely, 
the coefficient of h(j) i^A in equation f )58|) goes to 

1.4 3Q 2 

^^(2 + rT) (6i) 
9 

when /i — )• 00. 
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The remaining terms in equation (1581) have also finite limits when h — > oo. Finally, we 
get the following NLS equation valid for infinite depth and constant vorticity : 

iA T - " (1 + ^ )2 % = -^^'w'a + ^ U + en + en 2 + rf ) \AfA (62) 

P(2 + fi) 3 ?5 2 C 2 l + fi 1 1 2 C 2 V y 1 1 v ; 

This equation deals with A which is the value of the velocity potential for y = 0. Using 
equation ( 138|) . the NLS equation for the enveloppe of the wavetrain is : 



^(1 + ^) 2 u;fc 2 ^ 2 (2 + ^) 2 | / - -2 t=t3\ i , 2 

za r =— a« = — _ a a H 4 + 6i2 + 6i2 + U a a 63 

P(2 + fi) 3 ?€ 8 1 + Q 1 1 8 V 7 1 1 v ; 

We have left deliberately two nonlinear terms. The first term of the RHS comes from the 

coupling between the mean flow and the vorticity while the second can be obtained heuris- 

tically from the nonlinear dispersion relation. When Q vanishes, this coupling disappears 

and the heuristic method can be applied. Note that Baumstein^ and Li et al.— missed this 

coupling. 

If we use the heuristic method to obtain the NLS equation from the nonlinear dispersion 
relation that was found by Simmen & Saffman 4 , we should obtain only the second term of 
the RHS of (|63|) . So, we have shown that the heuristic method is not valid in presence of 
vorticity, even in infinite depth. 

Equation (|63|) is rewritten as follows : 

k 2 (2 + n) 3 ^ 8(i + n) 

The dispersive and nonlinear coefficients of equation (154"]) present two poles VL = —2 and 
Q = — 1 and two zeros Q = —1 and Q = —2/3 respectively. Nevertheless, we recall that 
Q > — 1 and consequently in infinite depth Q will never be equal to —1 or —2. 
For Q = —2/3, the nonlinear coefficient vanishes and the vor-NLS equation is reduced to a 
linear dispersive equation (the Schrodinger equation) 

3 

ia T - g^a« = (65) 
III. STABILITY ANALYSIS AND RESULTS 



ia T - U{1 + ^ 2 a^ = U + 10H + 8H 2 + 3rf) |a| 2 a (64) 

h? 2 + K s unn / 1 1 v ; 



The equation ( 1561) admits the following Stokes's wave solution 

a = a exp(— iMa^r) (66) 
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We consider the following infinitesimal perturbation of this solution 

a = a (l + 5 a ) exp[i(5 u) - Ma^r)] (67) 

Substituting this expression in equation ( 156|) and linearizing about the Stokes' wave solution, 
we obtain 

.dS a 96 u . 2, T d2s a . . T d 2 5 w 2 

- -q^T + 5 aMa + L— + iL— - 3Ma 5 a = (68) 

Separating the real and imaginary parts, the previous equation transforms into the following 

system 

d5 a d 2 5 w _ 
L 4 M ( 69 ) 

This is a system of linear differential equations with constant coefficients that admits the 
following solution 

5 a = A a exp[i(l£ - At)] 

(70) 

5 U = A w exp[z(/£ - At)] 
Substituting this solution in the system of equations (IBTJj) gives 

i\A a + l 2 LA w = 

(71) 

(2Ma 2 + l 2 L)A a — i\A w = 
The necessary and sufficient condition of non trivial solutions is : 

A 2 = £ 2 L(2Ma 2 + £ 2 L) (72) 

Discussion : When L{2Ma^ + £ 2 L) > there are two real solutions, the perturbation 
is bounded and the Stokes' wave solution is stable while when L(2Ma^ + l 2 L) < the 
perturbation is unbounded and the solution is unstable. Note that the latter condition 
implies that LM < 0. 

We set L = L\^2 and M = M\ujk 2 so that L\ and M\ are dimensionless functions of kh 
and Q only. The growth rate of instability is then 

Iuj 



l = 1 - 2 -\/-2M 1 L 1 k*a 2 -PL 2 (73) 



Its maximal value is obtained for I = y — j^-aok 2 and is : 

7max = M^ikao) 2 (74) 
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Figure 3. Stability diagram in the (0, fc/t)-plane. S : stable, U : unstable. 

The instability domain is plotted in figure [3] as a function of the parameters Q and kh. As 
soon as— 1/a < Q < —2/3 the waves become stable to modulational perturbations. Noting 
that Q is an increasing function of Q, it is easy to show that — 1/a < Q < —2/3 corresponds 
to — oo < Q < —2\J~^-. Hence, there is a value Q c = — ^\J~^§ of Q (depending on k) for which 
Q = —2/3. For a constant vorticity corresponding to Q < Q c waves are linearly stable. In 
particular, Stokes' waves of wavenumber k propagating on a linear shear current satisfying 
Q < —2\J~^- are stable to modulational instability whatever the value of the depth may be. 
There is a critical value kh cr a of the parameter kh, as shown in figure [31 above which 
instability prevails. For Q = (no vorticity) this threshold has the well known value 1.363. 
The critical value of this threshold is reached very near Q = 0. 

The linear stability of the Stokes wave solution is known to be controlled by the sign of the 
product LM of the coefficients of the vor-NLS equation floT)]) . Let us consider this product 
when kh — > oo 

8 (2 + fi) 3 1 ' 

The condition LM < corresponds to instability whereas LM > corresponds to stability. 
In the domain Q > —1, this product admits one simple root Q = —2/3. For this value of Q, 
LM changes sign and as a result there is an exchange of stability. Hence, in infinite depth 
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we can claim that there is no modulational instability when — 1 < Q < —2/3. 
In order to illustrate the restabilisation of the modulational instability we consider a mod- 
ulated wave packet that propagates initially without current in infinite depth and meets 
progressively a current with Q = —0.83 which corresponds to a stable regime. The results of 
the numerical simulations of the vor-NLS equation are shown in figures H] and [5j Temporal 
evolutions of the ratio A max (t)/A are plotted without (Q = 0) and with (Q = —0.83) shear 
current where A max (t) and Aq are the maximum amplitudes of the modulated wave train at 
time t and time t = 0, respectively. In figure H] the vorticity is initially set equal to zero. At 
t = 200 the value of Q is increased progressively up to —0.83 (that belongs to ] — 1, —2/3]) at 
t = 600 and remains equal to this value till the end of the numerical simulation. The carrier 
amplitude and carrier wavenumber are kA = 1/16 and k = 8 respectively. The perturbation 
amplitude is one tenth of the carrier amplitude and its wavenumber is Ak = 1 = 1. Hence, 
the criterion for the occurrence of a simple recurrence is satisfied. For Q = 0, one can ob- 
serve the Fermi-Pasta-Ulam recurrence phenomenon (FPU) which corresponds to a series of 
modulation-demodulation cycles. When the shear current is introduced the Benjamin-Feir 
instability is strongly reduced. In figure [5j the same numerical simulation is conducted, but 
the wave steepness of the carrier wave is now kAo = = 0.1083 and so the wavenumber 
21 corresponds to an unstable perturbation. In figure [5] is shown a double recurrence in the 
absence of shear current. The introduction of the vorticity modifies drastically this recur- 
rence. When Q reaches the value —0.83, the modulational instability is removed. Note in 
the presence of vorticity the increase of the amplitude of the envelope of the wave packet 
near t = 400. At this time Q does not yet belong to the stable interval ] — 1, —2/3]. 

A. Growth rate of instability 

The ratio of the maximum growth rate of instability given by equation (1741) to its value 
in the absence of shear is plotted in figure [6] as a function of Q for Q > —2/3 and several 
values of kh. In infinite depth, the presence of vorticity increases or decreases the maximum 
growth rate of modulational instability, 7 max , when Q > or — 2/3 < Q < 0, respectively. 
In finite depth and — 2/3 < Q < 0, the effect of vorticity is to reduce the maximum rate of 
growth whereas for O > we observe an increase and then a decrease. 

In figure [7] is shown the behavior of the normalized maximum growth rate as a function of 
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Figure 4. Temporal evolution of the normal- 
ized maximum amplitude of the envelope in 
the case of a simple recurrence for kh = oo : 
H = (solid line), U = -0.83 (dash-dotted 
line) 



Figure 5. Temporal evolution of the normal- 
ized maximum amplitude of the envelope in 
the case of a double recurrence for kh = oo : 
Ti = (solid line), H = -0.83 (dash-dotted 
line) 



kh for several values of Q. Herein, the normalization is different from that used in figure EJ 
Figure [7] correspond to values of Q larger than —2/3. For Q > —2/3, the critical value kh CT n 
associated to restabilisation is very close to 1.363 and corresponds to Q ~ 0. In figure [7] for 
Q > —2/3 the maximum growth rate of instability increases with kh greater than 1.363. 

In figure [8] is plotted the normalized rate of growth of modulational instability as a 
function of the perturbation wavenumber £ for several values of fl, within the framework of 
finite depth. Figure M corresponds to the case of infinite depth. 



B. Bandwidth instability 

In figure [10] is shown the ratio of the instability bandwidth A£ to its value in the absence 
of shear current A£q = A£(Q = 0) as a function of Q for several values of kh. From equation 

k 2 ao. One can observe an increase of the band of 
instability followed by a decrease when f2 increases, except when depth becomes infinite. 

In table 1 is presented a comparison of our results with those of Oikawa et al. (1987) 
in the case of two dimensional flows for two values of kh and several values of the Froude 
number. Note that the Froude number, F, they used is exactly Q. This comparison shows 
a quite good agreement between Oikawa et al. and present results. 
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Figure 6. Normalized maximum growth rate as a function of ft for kh = 1.40 (solid line), kh = 1.70 
(dashed line) and kh = oo (dash-dotted line). 7o max is the maximum growth rate in the absence of 
shear current 
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Figure 7. Normalized maximum growth rate as a function of kh for Q = (solid line), Q = —0.50 
(dashed line) and Q = 3.0 (dash-dotted line). 7o max is the maximum growth rate when kh = oo 
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Figure 8. Normalized growth rate func- 
tion of the perturbation wavenumber t for 
kh = 2.0 and Tl = -0.5 (solid line), H = 0.0 
(dashed line), Q = 0.5 (dot-dashed line) 



Figure 9. Normalized growth rate func- 
tion of the perturbation wavenumber I for 
kh = oo and = —0.5 (solid line), f2 = 0.0 
(dashed line), Q = 0.5 (dot-dashed line) 





F = 0.0 


F = 0.25 
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F = 2.0 
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1.0/1.00 


1.2/1.31 


6.0/5.96 
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Table I. Comparison with results of Oikawa et al. (1987) : F is the Froude number. The first value 
is estimated from their figures whereas the second one corresponds to our computations with the 
vor-NLS equation 



C. Benjamin- Feir index in the presence of vorticity : Application to rogue 
waves 

Within the framework of random waves Janssen^ introduced the concept of the Benjamin- 
Feir Index (BFI) which is the ratio of the mean square slope to the normalized width of the 
spectrum. When this parameter is larger than one, the random wave field is modulationally 
unstable, otherwise it is modulationally stable. From the NLS equation Onorato et a/.— 
define the BFI as follows 



(76) 
17 



BFI 



aok 
Ak/k 



M 1 



7 
6 - 
5 - 



o . 
n3 4 

CM 

o 

^ 3 h 




Figure 10. Normalized instability bandwidth as a function of f2 for kh = 1.5 (solid line), kh = 1.8 
(dashed line), kh = oo (dot-dashed line) 

where Ak represents a typical spectral bandwidth. 
In infinite depth the BFI without shear current is 

4a k 



BFIn 



Ak/k 



(77) 



Hence, the normalized BFI writes 



BFI 



M 1 



(78) 



BFI, 

The coefficients Mi and L\ depend on the depth and vorticity. Onorato et a/.— considered 
the effect of the depth on the BFI. Herein, besides depth effect a particular attention is paid 
on the influence of the vorticity on the BFI. In order to measure the vorticity effect on the 
BFI, the ratio of the BFI in the presence of vorticity to its value in the absence of vorticity 
in infinite depth is plotted in figures [TTJ and [12j For fixed value of Q the BFI increases 
with depth. Our results for Q = are in full agreement with those of Onorato et al. 20 (the 
solid line in figure ITT1) . Furthermore, it is shown for Q > and sufficiently deep water that 
the BFI increases with the magnitude of the vorticity. Therefore, we may expect that the 
number of rogue waves increases in the presence of shear currents co-flowing with the waves. 
For Q < the presence of vorticity decreases the BFI. For a more complete information 



18 



1 23456789 10 
kh 

Figure 11. Normalized Benjamin Feir Index 
as a function of kh for several values of f2 : 
O = 0.0 (solid line), Q = 1.0 (dashed line), 
12 = 2.0 (dot-dashed line) 
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kh 

Figure 12. Normalized Benjamin Feir Index 
as a function of kh for several values of f2 : 
17 = 0.0 (solid line), Q = —0.3 (dashed line, 
12 = —0.6 (dot-dashed line) 



about rogue waves, one may consult Kharif, Pelinovsky and Slunyaev (2009)^. 



IV. CONCLUSION 

Using the method of multiple scales, a ID nonlinear Schrodinger equation has been 
derived in the presence of a shear current of non zero constant vorticity in arbitrary depth. 
When the vorticity vanishes, the classical NLS equation is found. A stability analysis has 
been developed and the results agree with those of Oikawa et al. (1987) in the case of 
ID NLS equation. We found that linear shear current may modify significantly the linear 
stability properties of weakly nonlinear Stokes waves. 

We have shown the importance of the coupling between the mean flow induced by the 
modulation and the vorticity. This coupling has been missed (or not emphasized) by previous 
authors. 

Furthermore we have shown that the Benjamin- Feir instability can vanish in the presence 
of positive vorticity (Q < 0) for any depth. 
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